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Abstract 

In the virtual presence of a heavy quark t, the interactions of a CP-odd scalar 
boson A, with mass Ma <C 2M^, with gluons and light quarks can be described by an 
effective Lagrangian. We analytically derive the coefficient functions of the respec- 
tive physical operators to three loops in quantum chromodynamics (QCD), adopting 
the modified minimal-subtraction (MS) scheme of dimensional regularization. Spe- 
cial attention is paid to the proper treatment of the 75 matrix and the Levi-Civita e 
tensor in D dimensions. In the case of the effective gg A coupling, we find agreement 
with an all-order prediction based on a low-energy theorem in connection with the 
Adler-Bardeen non-renormalization theorem. This effective Lagrangian allows us to 
analytically evaluate the next-to-leading QCD correction to the A — > gg partial de- 
cay width by considering massless diagrams. For Ma = 100 GeV, the resulting cor- 
rection factor reads l + (221/12)ai 5) (M A )/7r+165.9 (cSP(M A )/n) « 1+0.68+0.23. 
We compare this result with predictions based on various scale-optimization meth- 
ods. 

PACS numbers: 11.15.Me, 12.38.Bx, 13.30.Eg, 14.80.Cp 

*Permanent address: Institute for Nuclear Research, Russian Academy of Sciences, 60th October 
Anniversary Prospect 7a, Moscow 117312, Russia. 

T Permanent address: Max-Planck-Institut fiir Physik (Werner-Heisenberg- Institut), Fohringer Ring 6, 
80805 Munich, Germany. 



1 Introduction 



Despite the tremendously successful consolidation of the standard model (SM) of ele- 
mentary particle physics by experimental precision tests during the past few years, the 
structure of the Higgs sector has essentially remained unexplored, so that there is still 
plenty of room for extensions. A phenomenologically interesting extension of the SM 
Higgs sector that keeps the electroweak p parameter [|IJ at unity in the Born approxima- 
tion, is obtained by adding a second complex isospin-doublet scalar field with opposite 
hypercharge. This leads to the two-Higgs-doublet model (2HDM). After the three mass- 
less Goldstone bosons which emerge via the electroweak symmetry breaking are eaten up 
to become the longitudinal degrees of freedom of the W ± and Z bosons, there remain 
five physical Higgs scalars: the neutral CP-even h and H bosons, the neutral CP-odd A 
boson, and the charged if ± -boson pair. The Higgs sector of the minimal supersymmetric 
extension of the SM (MSSM) consists of such a 2HDM. At tree level, the MSSM Higgs 
sector has two free parameters, which are usually taken to be the mass Ma of the A 
boson and the ratio tan/3 = v^/vi of the vacuum expectation values of the two Higgs 
doublets. For large values of tan /3, the top Yukawa couplings of the neutral Higgs bosons, 
$ = h,H, A, are suppressed compared to the bottom ones. 

The search for Higgs bosons and the study of their properties are among the prime 
objectives of the Large Hadron Collider (LHC), a proton-proton colliding-beam facility 
with centre-of-mass energy yfs = 14 GeV, which is presently under construction at CERN. 
The dominant production mechanisms for the neutral Higgs bosons at the LHC will be 
gluon fusion, gg — > $ 0, and bb& associated production, gg,qq — > bb& ||, which is, 
however, only relevant for large tan/?. The loop-induced gg& couplings @] are mainly 
mediated by virtual top quarks, unless tan /3 is very large, in which case the bottom- 
quark loops take over. The ggh and ggH couplings also receive contributions from squark 
loops, which are, however, insignificant for squark masses in excess of about 500 GeV 
||. In the case of the g g A coupling, such contributions do not occur at one loop because 
the A boson has no tree-level couplings to squarks. For small tan j3, the inclusive cross 
sections of pp — > $ + X via gluon fusion are significantly increased, by typically 50-70% 
under LHC conditions, by including their leading QCD corrections, which involve two- 
loop contributions || 0]. Thus, the theoretical predictions for these observables cannot yet 
be considered to be well under control, and it is desirable to compute the next-to-leading 
QCD corrections at three loops, since there is no reason to expect them to be negligible. 
Recently, a first step in this direction has been taken by considering the resummation of 
soft-gluon radiation in pp — > $ + X, assuming tan/3 to be small ||. 

Considering the enormous complexity of the exact expressions for the leading QCD 
corrections 0, it becomes apparent that, with presently available technology, the next- 
to-leading corrections are only tractable in limiting cases. For instance, in the large-tan /3 
limit, where the gg<& couplings are chiefly generated by bottom-quark loops, one can 
neglect the bottom-quark mass against the Higgs-boson mass, keeping the bottom Yukawa 
coupling finite. In this way, one resorts to massless QCD. On the other hand, if tan/3 
is close to unity, which was assumed in Ref. ||, the top-quark loops play the dominant 
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role, and simplifications occur if the mass hierarchy M$ <C 2m t is satisfied. Then, it 
is useful to construct an effective Lagrangian by integrating out the top quark. This 
effective Lagrangian is a linear combination of local composite operators of mass dimension 
four, which act in QCD with five massless quark flavours, while all dependence on mt is 
contained in their coefficient functions. Once the coefficient functions are known, it is 
sufficient to deal with massless Feynman diagrams. The effective Lagrangian describing 
the interactions of the SM Higgs boson H with gluons and light quarks was elaborated at 



two loops in Ref. |J and extended to three loops in Refs. (10], |TT| . As an application, the 
0{a s ) H and 0(a 2 ) [|Kj corrections to the H — > gg partial decay width were calculated 
from this Lagrangian. These results can be readily adapted to the h and H bosons of the 
2HDM with small tan (3 by accordingly adjusting the top Yukawa coupling, which appears 
as an overall factor. 



In this paper, we extend the three-loop analysis of Refs. |ll]] to include the A boson 
of the 2HDM. As in Ref. we work in the limit tan (3 ~ 1, so that we may treat bottom as 
a massless quark flavour with vanishing Yukawa coupling, on the same footing as up, down, 
strange, and charm. Specifically, we construct a heavy-top-quark effective Lagrangian for 
the QCD interactions of the A boson and derive from it an analytic result for the 0(a 2 ) 
correction to the A —>■ gg partial decay width appropriate for Ma <C 2m t . We recover 
the corresponding 0(a s ) result originally found in Refs. |7], |12|] and also discussed in 
Ref. fl3||. The 0(a s ) correction for arbitrary values of the A-boson and quark masses was 
presented in Ref. as a two-fold parameter integral. Furthermore, it was shown that the 
leading high-m t term of this correction may also be obtained from massless five-flavour 
QCD endowed with a heavy-top-quark effective ggA coupling J7], |TJ]]. A central ingredient 
for this check was the observation that the effective ggA coupling does not receive QCD 
corrections, at least at 0(a s ). This fact was interpreted [J7], [TJJ as being a consequence 
of the Adler-Bardeen theorem fl5|, which states that the anomaly of the axial- vector 



current |16| is not renormalized in QCD. This theorem is strictly proven to all orders in 



a s for the abelian case [15], and strong arguments suggest that it also holds true for the 



nonabelian case ||17|| . In this paper, we verify by an explicit diagrammatic calculation 
that the 0(a s ) and 0(a 2 ) corrections to the coefficient function of the operator G a ^ v G a ^ v , 
which generates the A-boson effective couplings to gluons vanish.^] We also present the 
leading-order coefficient function of the physical operator pertaining to the effective qqA 
interaction, where q is a light quark. We thus provide the tools which are necessary 
to reduce the calculation of the next-to-leading QCD correction to the cross section of 
pp — > A + X to a standard problem in massless five-flavour QCD. 

In our analysis, we consistently neglect the Yukawa couplings of the light quarks to 
the A boson. In other words, if all quark masses, except for mt, are nullified, the hadronic 
decay width of the A boson is entirely due to A — > gg and the associated higher-order 
processes under consideration here. Through three loops, the contributing final states are 

9999, 99W, qqq'q', 999, gqq, 99, and qq. 

This paper is organized as follows. In Section 0, we establish the heavy-top-quark 



1 Strictly speaking, this statement is only true as long as we ignore the axial-anomaly equation, as will 
become apparent in Section 0. 
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effective Lagrangian for the QCD interactions of the A boson to three loops. In Section 
we compute from this Lagrangian the 0(a 2 s ) correction to the partial width of the decay 
A — > gg and compare it with predictions based on various scale-optimization methods. 
Section f| contains a discussion of our results together with some remarks on the connec- 
tion between the effective ggA coupling, the axial-anomaly equation, and the low-energy 
theorem. 



2 Effective Lagrangian 

We start by setting up the theoretical framework for our analysis. As usual, we em- 
ploy dimensional regularization in D = 4 — 2e space-time dimensions and introduce a 



't Hooft mass, to keep the coupling constants dimensionless [18, 19]. We perform the 
renormalization according to the modified minimal-subtraction |2T| (MS) scheme. 
For the sake of generality, we take the QCD gauge group to be SU(iV c ), with N c ar- 
bitrary. The adjoint representation has dimension Na = N% — 1. The colour factors 
corresponding to the Casimir operators of the fundamental and adjoint representations 
are Cp = (N% — 1)/(2N C ) and Ca = N c , respectively. For the numerical evaluation, we 
set N c = 3. The trace normalization of the fundamental representation is T = 1/2. As 
an idealized situation, we consider QCD with m — rif — 1 light quark flavours and one 
heavy flavour t, in the sense that 2m qi -C Ma -C 2mt. We wish to construct an effective 
nj-flavour theory by integrating out the t quark. We mark the quantities of the effective 
theory by a prime. Bare quantities carry the superscript "0". As already mentioned in the 
Introduction, we consider a 2HDM with tan/3 = 1, so that the quark Yukawa couplings 
and masses are related by a flavour-independent proportionality factor. 

The starting point of our consideration is the bare Yukawa Lagrangian for the inter- 
actions of the A boson with the quarks in the full nj-flavour theory, 

C = -^(t + , (1) 

where v = 2- 1 / 4 G F 1/2 , with G F being Fermi's constant. Taking the limit m" — > oo and 
keeping only leading terms, Eq. ([]]) may be written as a linear combination of pseudoscalar 
composite operators, O^, with mass dimension four acting in the effective ni- flavour theory. 
The resulting bare Lagrangian reads 

4° 

c eS = -—(c Q 1 d' 1 + cld' 2 + ..) ) (2) 

where 

qi qOI ,a qQI ,afiu 

ni 

6' 2 = d,J ^, Jf^E ^75?°', (3) 

i=i 
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Cf are coefficient functions, which depend on the bare parameters of the full theory and 
carry all dependence, and the ellipsis stands for terms involving unphysical operators, 
which do not contribute to physical observables. Here, G°L, = d^G^ — d v G" + gsf^G^Gl 
is the colour-field-strength tensor and G apu = e pvpa G a pcF is its dual; G a (a = 1, . . . , N A ) 
are the gluon fields, g s = y/4:ira s is the QCD gauge coupling, and f abc are the structure 
constants of the SU(iV c ) algebra. We do not display the colour indices of the quark fields. 
We mark the operators and coefficient functions with a tilde in order to avoid confusion 



with our previous notation for the scalar case pi, 11 



The Levi-Civita tensor e p,l,pa is unavoidably a four-dimensional object and should be 
taken outside the R operation. Thus, we rewrite Eq. (Q) as 0[ 



22, 23 



AivpoQi where 

c ^l, pupal WUC1C 



U l,pvpa — (jr [pv (j pa] ' 



n, 



o' 



2, pupa 



5i™ 



^0/ 



■1 



[pQi ivipi^n 



0/ 



(4) 



are antisymmetrized in their four D-dimensional Lorentz indices. Furthermore, we substi- 
tute 75 = (z/4!)e Atlyp °7[ Ai 7 ; ,7p7 (J ] in Eq. (JI]). We then carry out the D-dimensional calcula- 
tions with e pupcr peeled off from the expressions. In the very end, after the renormalization 
is performed and the physical limit e — > is taken, we contract the expressions with e^ upa 
to obtain the final results. 

Prior to describing the actual calculation of C®, let us discuss how Eq. (|2|) is renormal- 
ized. Since we are only interested in pure QCD corrections, we may substitute A = A 
and v° = v in Eqs. ([!]) and (§). Denoting the renormalized counterparts of Cf and 0[ by 
Ci and [O'j], the renormalized version of Eq. (|2]) takes the form 



= -2 l '^GfA (Ct [d' 1 ]+C 2 [d' 2 



+ 



(5) 



5 /,Al is renormal- 
ized multiplicatively in the same way as the colour- singlet axial- vector current Jg itself, 
while G f l^G '' apu mixes under renormalization [ 



where the ellipsis again represents unphysical terms. The divergence d p J ; 

det axial- vector 
24j. Specifically, we have 



o;] = z' n d[ + z[ 2 d' 2 , 

2 = Z 22 2 , 



(6) 



where Z 22 



Zy^Z^ is the product of the standard ultraviolet (UV) renormalization 



constant Zf^ s of the singlet axial current in the MS scheme and the finite renormalization 
constant Zf. The latter is in 
relation of the axial anomaly, 



constant Zf. The latter is introduced to restore the one-loop character of the operator 



O' 



a 



7T 



(aO Tni 
4 



(7) 
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Figure 1: Typical Feynman diagrams contributing to the coefficients C® in Eq. (|2]). 
Looped, bold-faced, and dashed lines represent gluons, t quarks, and A bosons, respec- 
tively. 



which is valid for Pauli-Villars regularization ||23|| . For the reader's convenience, we list 
here the various renormalization constants to the order necessary for our purposes. They 



read 23 



Z' u = l + 



a 



(//) 1 



%12 
MS 



X 

a 



7r e 
1 /121 



Y 2 c a + \tu i] + 



fry 



r 2 

M44° A 



11 
18 



1, 
9 



CaTth + ^T 2 ^) + 



17 C\ + \c F Tm + ^-CaTti, 



18 



24 



7T 



-(3CV) 



i. 



Note that Z' n = {z'^j is the square of the coupling renormalization constant, Z' g . As will 
become apparent later, we only need the leading term of Z' 22 . The relations between the 
bare and the renormalized coefficient functions are accordingly given by 



c 2 



1 ^ 
1 

Z'22 



(9) 



We now turn to the computation of the bare coefficient functions and C 2 in Eq. (|2]). 
We are thus led to consider irreducible vertex diagrams which connect one A boson to the 
respective states of gluons and light-quark pairs via one or more t-quark loops, whereby 
all external particles are taken to be on their mass shells. Typical examples are depicted 
in Fig. [I]. There are three independent ways to obtain C°, namely from the sets of gg A 
three-point, gggA four-point, or ggggA five-point diagrams. At the three- loop level, these 
sets contain 657, 7362, and 95004 diagrams, respectively. We choose to work out the first 
option in the covariant gauge with arbitrary gauge parameter, so that the gauge-parameter 
independence of the final result yields a nontrivial check. Another independent check is 
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then provided by elaborating the second option in the 't Hooft-Feynman gauge keeping 
only one external momentum different from zero. In order to cope with the enormous 
complexity of the problem at hand, we make successive use of powerful symbolic ma- 
nipulation programs. Specifically, we generate and evaluate the contributing diagrams 
with the packages QGRAF |25| and MATAD |2§, which is written in FORM fSTJ, respec- 
tively. The cancellation of the UV singularities, the gauge-parameter independence, and 
the renormalizat ion-group (RG) invariance serve as strong checks for our calculation. 

Let us denote the sum of all relevant ggA diagrams by V^° af3 (qi, q 2 ), where qi and q 2 
are the incoming four-momenta of the two on-shell gluons with polarization four- vectors 
and respectively. According to Eq. (f|), V^^ aj3 is by construction totally antisymmetric 
in the indices /z, v, p, and a. In order to compute C°, we need to expand V gg A up to 
terms linear in q\ and q 2 - There is just one possible structure, namely 



(10) 



so that the coefficient C gg A may be conveniently extracted by noting that P gg A = —{D — 
2)(D - 3)(gi ■ g 2 ) 2 /24. The final formula for C° reads 



7P 7P 



(D-2)(D-3)( qi -q 2 ) 2 Z, 



A3 



PggA^vpailli ^)^ gg A,aM x ' ^ L 



:n) 



The normalization of Eq. may be understood as follows. On the one hand, we need 
to renormalize the mass and the pseudoscalar current of the t quark in the Lagrangian ([!]) 
of the ^/-flavour theory. In this way, the renormalization constants Z, 
23] enter. They are defined as 



28 



-^Ms> ari d Z\ 



[i l5 t] = Z^Zif^t 



;i2) 



and read 23, 28 



7 P 



(n/ V) 




07 k 

— C F C A + —C F Trif 
192 48 S 



1 + 



+ 



a 



(»/) 



(/i) 1 



Zl = 1 + 



a 



7T e 

~64° F 

(B/) (aO 



TV 



79 11 

C F C A C F Trif 

192 48 } 



-2C F ) + 



o 



( n f) 



7T 



7T 



32 



32 



-C F Tn 



— ( — Cp + —CpCa 



32 



32 



-CpTn 



^CfCa + ^C F Tn f 



(13) 
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The finite renormalization constant Z\ is needed in addition to the usual UV renormaliza- 
tion constant of the MS scheme, Z^, to effectively restore the anticommutativity of the 
75 matrix [^3, 29|. On the other hand, we need to express the bare couplings and fields 



appearing in the Lagrangian (g) of the ^-flavour theory in terms of their counterparts in 
the rij-flavour theory. The appropriate relations, 

9? = (° 9 9°s, <f = Vfe G ^ = ( 14 ) 

involve the decoupling constants C° > an d C3 • m the case °f the ggA amplitude, only £° 
occurs. For our purposes, we need £3 through 0(e 2 a s ) and 0(ea 2 s ). The corresponding 
expression may be extracted from Eq. (B.2) of Ref. ||11|| , where the renormalized version, 



( 3 , is listed through 0{a*) in the covariant gauge. Finally, a factor of 8 stems from the 
Feynman rule for the two-gluon piece of 0[. In Eq. (]TTj) , we may take the limit (ft, tfe ~~ > ; 
which reduces the problem of finding C° to the solution of massive vacuum integrals 
After renormalization according to Eq. (|9|), we find 



71 



1 1 Q ' w ' )(/x) I ( ^M )' 



IT \ 7T 



(15) 



i.e. the correction terms of C(a s ) and 0(a 2 s ) indeed vanish, as was suggested in Ref. 0, [14 



on the basis of the nonabelian variant [17] of the Adler-Bardeen theorem [^5| , which 



predicts that Ci does not receive QCD corrections. Notice that this is only true if C% is 
expanded in a^ n ^(/i). 

As an independent check, we may also extract C° from the sum VggJ% a p Jq) of the 
gggA diagrams, where a, (3, and 7 are the Lorentz indices of the gluon polarization four- 
vectors. Notice that it is sufficient to keep one external four-momentum, q, different from 
zero, since the three-gluon piece of 0[ contains just one derivative. Again, there is only 
one possible structure linear in q, namely 

VgggA,al3'y{ ( l) = ^ gggAPgggA^^W) i ^ggg^afijil) = Q M '9a9p9°\ (16) 

so that the coefficient C ggg A may be easily extracted by using Pg gg A — (D — 1) (D — 2) (D — 
3)g 2 /24. In contrast to the two-gluon piece of 0[, we now have to include decoupling 
constants for three gluon fields and one gauge coupling. Again, the Feynman rule for the 
three-gluon piece of 0[ involves a factor of 8. Thus, the final formula for C® is given by 

(jO _ 2 Zj MS Zj 5 pa(37 ( n\V^ vpcr (a) (17) 

1 (D- l)(D-2)(D-3)g 2 Z m C 5 (C 3 °) 3/2 999Aq/37 

Taking the limit q — > and and working in the 't Hooft-Feynman gauge, we recover from 
Eq. ( |P7D our previous result for (7°. 

Finally, we turn to (7°, which is generated by vertex diagrams that couple the A 
boson to a pair of light quarks and involve a virtual t quark. This requires at least 
two loops. Specifically, there are 2 two-loop and 63 three-loop diagrams of this type. 
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Calling the resulting qiqiA amplitude V^ q u / a (q) , where the argument q is the incoming 
four-momentum of the A boson, we may extract C 2 as 



a yV 7P 

2 p-l)(£>-2)p-3)g 2 Z m C 2 ° L^ 7 ^^ W 



In order to treat the A —>■ gg decay at three loops, it is sufficient to know the leading 
term of C^- Moreover, due to Eq. (||), the computation of the next-to-leading term 
of C2 would require the knowledge of Z' 12 to 0(a 2 ), which is not yet available. After 
renormalization according to Eq. ([|), we find 

V t / \16 8 mUtt))' K ' 

where m t (/i) is the MS mass of the t quark. 



3 A — ► gg decay 

Having established the high-m 4 effective Lagrangian ([]) controlling the QCD interactions 
of the A boson, we are now in a position to evaluate from it the 0(a 2 ) correction to 
the A — > gg decay width. To this end, we need to compute the absorptive part of the 
A-boson self-energy, at q 2 = M|, induced by [0[] and [0' 2 ] to sufficiently high order in 
the nj-flavour theory. 

In D dimensions, the correlator function, at four-momentum q, of two bare operators 
of the type defined in Eq. (f|) has the Lorentz decomposition 

(oU^r VpV ) = K^^fJg^'g^ + ^WmVffZaigQ (20) 

where II^ and II^,- are functions of q 2 . We may extract 11° ^ and II^,- by totally 
contracting Eq. (|20|) with the projectors 

24 [q 9 [fl 9 U 9p 9 a ] ~ 9u 9p 9 a ] 



Pi, pupa (l) 



(q 2 ) 2 (D- l)(D-2)(D-3)(D-4) 

, 96 (-q 2 g [ £gvg p p'g$ + Dg^' g{ 9$ , 

2,pupa W (<?2)2 p_l)(£)_2)(£)-3)(D-4) ' 



respectively. Notice that Pi and P 2 develop 1/e poles in the physical limit e — > 0. This 
may be understood by observing that the two terms on the right-hand side of Eq. fl2D| ) are 
actually linearly dependent in four dimensions. In practice, the appearance of 1/e poles 
in Eq. ([H]) does not create a problem, since we are only interested in the absorptive parts 
of the correlators in Eq. (|20|), so that it suffices to extract the pole parts of the relevant 
diagrams. It turns out that = 0, so that, after performing renormalization, taking 
the physical limit e — > 0, and contracting with the Levi-Civita tensors, we have 

31 fO'l) = -6g 2 n^(g 2 ), (22) 
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Figure 2: Typical Feynman diagrams contributing to the correlator (0[0[). Looped, 
solid, and dashed lines represent gluons, light quarks, and A bosons, respectively. Solid 
circles represent insertions of 0[. 



where n 2 ,ij is the renormalized version of 11° j,-. In fact, it can be shown that n^,- van- 
ishes on kinematical grounds. As is well known, 0[ can be written as the divergence 



of the so-called Chern-Simons current, K' ,(Jl 



AG° u '' a d p G' 



Or, a 



+ 



i.e. 0[ = d^K''^, which is an exact identity. This implies that 
0'i,pu P a = dlfj.K'vpa]- Thus, the correlator in Eq. ( p0|) is represented by just one term 
proportional to q^q^ 1 (K' ■<K'>' /pa ^), whence it follows that 11°^ = 0. 

Due to Eq. (§), all three correlators (d[d[), (d[6 2 ), and (0 2 6 2 ) contribute to ([0£][0£]), 
the absorptive part of which we wish to calculate through 0(a 2 ). At the three-loop level, 
these three correlators receive contributions from 403, 28, and 33 massless diagrams, re- 
spectively. Typical examples pertaining to (0^0^) are depicted in Fig. |^. We generate 
and evaluate the contributing diagrams with the packages QGRAF |25| and MINCER 
31"| , which is written in FORM [27|. We work in the covariant gauge with arbitrary 
gauge parameter. The cancellation of the latter in the final results serves as a welcome 
check. 

Our results for the absorptive parts of the renormalized correlators read 
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-N A < 1 
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(ni) 
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7T 
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0' 2 0' 2 



3 q 2 



+ 



7T 



1135, /i 2 121 2 /i 
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55 

C(2) - -C(3) 
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+ 3C(3)-21n^- 
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-^-n* + — n\ In — 



12 



Q- 



(g 2 ) 2 /JV A Tn, 



7T 



7T 



V 4 



(g 2 ) 2 (aW(ji)Y N A T 2 nf 



71 



71 



16 



(23) 



where ( is Riemann's zeta function, with values £(2) = 7r 2 /6 and £(3) ~ 1.202 057, and we 
have put N c = 3 in the second expression on the right-hand side. Notice that Im ([0 2 ][0 2 ]) 
starts at 0(a 2 ). This may be understood by observing that the A — > q^ decay width is 
helicity suppressed and quenched if the quark qi is taken to be massless. Thus, in order 
for a diagram to contribute it must have a cut which only involves gluons. Such diagrams 
first appear at three loops. Actually, the results in Eq. fl23|) are not mutually independent. 
In fact, Eq. © allows us to derive Im and Im ([0 2 ][0 2 ]) from Im ([Oi][Oi]) and 

thus provides a welcome check for Eq. The 0(a s ) term of Im ([O^O'J) in Eq. 



is in agreement with Ref. [12 



All ingredients which enter the calculation of the A — > gg decay width through 0{oQ] 
are now available. From Eq. (H), we derive the general expression 



m - gg) 



V2Gj 



C 2 Im ([Oil \0[]) + 2C 1 C 2 Im (\0[\ \6' 2 



+ Cflm 



O' O' 



(24) 

where it is understood that the correlators are to be evaluated at q 2 = M\. The last term 
contained within the parenthesis of Eq. ( f24"D contributes in 0{oP s ) and is only included for 
completeness. Inserting Eqs. fllED and (|2"3f) in Eq. (|24"D and putting N c = 3, rii = 5, and 
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jj, = M A , we finally obtain 



(5) 



71 



/237311 



x 



V 864 



(5) 



G F M\ ( 



- ^C(2) 
24 si y 

(M A ) N 




1 + 18.417 a ^ MA U 



x 171.544 - 5 In 



(25) 



The correction in Eq. (|2"5| ) agrees with the result originally found in Refs. [|7], [E|. 

If we assume that «( 5) (M A ) = 0.116, which follows from a^(M z ) = 0.118 for M A = 
100 GeV, and take the t-quark pole mass to be M t = 175.6 GeV, then the correction 
factor corresponding to the square bracket in Eq. (p5|) has the value 1 + 0.68 + 0.23 = 1.91, 
i.e. the three-loop term amounts to 33% of the two-loop term. This is somewhat larger 
than the corresponding correction factor for a SM Higgs boson with mass Mh = 100 GeV, 
which was found to be 1 + 0.66 + 0.21 = 1.87 ||10|| . For our choice of input parameters, we 
obtain from Eq. ( pop the QCD-corrected prediction T(A — > gg) = 426 keV. 

Similarly to the H —>■ gg case ||10|| , Eq. (^) may be RG improved by resumming the 
logarithms of the type ln(M t 2 /M^). This leads to 

2 



m - gg) 



a 



(• r >) 



G F M\ ( 

16ttv^ I * 



1 + 



4363 ctf(M A 



276 



7T 



60 af\M t ) 
23 



7T 



+ 



/239471 529 



C(2) 



V 864 24 

1800 af\M A ) af\Mi 
^29~ 



7T 



7T 




G F M\ (af\M A 
16ttV2 V 

+ 174.044 



+ 0.903 




1 + 15.808 s { — 



+ 2.609 



a: 



(6) 



{M t 



71 



71 



3.403 



af\M A ) af\M t 



7Y 



71 



(26) 



we 



Since, at present, the experimental lower bound on M A in the MSSM is 24.3 GeV [32 
have ln(M 2 /M|) < 4, so that the numerical effect of the RG improvement is negligible in 
practical applications. 

It is interesting to compare the exact value of the 0{a 2 s ) correction in Eq. (|25|) with 
the estimates one may derive from the knowledge of the 0(a s ) correction through the 



12 



application of well-known scale-optimization procedures, based on the fastest apparent 
convergence (FAC) ||33|| , the principle of minimal sensitivity (PMS) [0, and the proposal 
by Brodsky, Lepage, and Mackenzie (BLM) ]35| to resum the leading light-quark contri- 
bution to the renormalization of the strong coupling constant. These procedures lead to 
the generic expression 



71 



71 



G F M\ ( af\M A ) 



7T 



7[ 



(27) 



where K\ — k\ + niK%, with k\ = 97/4 and K\ = —7/6, is the coefficient of the 0(a s 
correction in Eq. (p5|) . The FAC, PMS, and BLM expressions for £, K\, and K 2 read 



-FAC 



lnf 

ine PMS 



Ki 



FAC 



0. 



2(3 I 2 3(3 0/ 



ln£ 



BLM 



-Kl, 



BLM 



33 

h + —Ki, 



PMS 



2A 
'Wo' 



K. 



PMS 



i? 2 BLM = Kl 



3,1/, 19\ 9 / , 363 
_« l7l a + _ f 3fcl + _j n ; - - llUi + — K! + 17 



:28i 



respectively, where /3o = 11/4 — and (3\ = 51/8 — 19ri;/24 are the first two coefficients 
of the Callan-Symanzik beta function of QCD. The numerical results for rti = 5 are 
summarized in Table 0. The values of K2 should be compared with the true coefficient 
K2 of the 0(a 2 s ) correction in Eq. (|25|). For completeness, we also list in Table [I] the 
corresponding results for the H — > gg decay width. In this case, one has K 2 ~ 156.808 — 
5.708 hi(Mf/Mfl) ]TD|. In both cases, all three scale-optimization prescriptions correctly 
predict the sign and the order of magnitude of K 2 - Furthermore, the three K 2 values for 
the H — > gg decay width are indeed smaller than the respective values for the A — > gg 
case. Similarly to Ref. 



36 , the FAC and PMS results almost coincide. 



4 Discussion and conclusions 

In this paper, we studied the interactions of a neutral CP-odd scalar boson A with gluons 
and Hi light quarks in the presence of a heavy quark t, with mass m t 3> Ma/2, through 
three loops in QCD. For simplicity, we assumed that the Yukawa couplings of the light 
quarks may be neglected, which is a useful approximation in the 2HDM with tan/3 w 
1, where the Yukawa couplings and the masses of the quarks are related by a flavour- 
independent proportionality factor. Starting from the Yukawa Lagrangian ([D) embedded 
in full QCD, we integrated out the t quark to obtain the corresponding effective Lagrangian 
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Table 1: Numerical evaluation of Eq. ( |28"D with m — 5 for the A — > gg and H ^ gg decay 
widths. 







H^gg 




Ki 


K 2 




Kt 


K 2 


FAC 


0.091 





277.601 


0.097 





263.346 


PMS 


0.081 


-0.841 


278.131 


0.087 


-0.841 


263.876 


BLM 


0.174 


5 


252.547 


0.174 


4.5 


242.484 



of the ni- flavour theory, the renormalized version of which is given by Eq. (|5|). The 
operators in Eq. (|5|) comprise only light fields, while all residual dependence on the t 
quark is contained in their coefficient functions. We diagrammatically evaluated the 
coefficients C\ and C 2 of the physical operators [0[] and [0' 2 ] through 0(oP s ) and 0(a 2 ), 
respectively. This is consistent, since [0' 2 ] is of 0(a s ) relative to [0[] as may be seen 
from Eq. (Uh. We worked in the MS renormalization scheme [B(J with the convention 
that Eq. (|7D should be exact f23|. Our results for C\ and C 2 are given in Eqs. (p~5|) and 
(P^D, respectively. In particular, we found that the 0(a s ) and 0(a 2 ) corrections to G\ 
exactly vanish if the latter is expressed in terms of a^ 1 ' (/i) . We thus verified, by explicit 
calculation through three loops, the all-order prediction that C\ does not receive any QCD 
corrections, which follows via a low-energy theorem |7], [14j from the nonabelian version 
1 17f| of the Adler-Bardeen non-renormalization theorem Jl5]. 

At this point, we should emphasize that, due to Eq. (0), the distinction between the 
renormalized operators [0[] and [0' 2 ] does not have a deep physical meaning. Actually, 
from Eqs. (|5|) and (|7|) it follows that the physical coupling of the A boson to [0[] is 
proportional to [C\ + a ( -™ t \ii)TniC 2 / (An)}. Thus, the low-energy theorem gets violated 
starting from 0{a\) unless we stick to the MS definitions of [0[] and [0' 2 ] |23[ . 



The effective Lagrangian (^j) allows us to evaluate physical observables related to the 
interactions of the A boson with gluons and light quarks to higher orders in QCD by just 
considering massless diagrams in the nj-flavour theory. As an application, we evaluated 
the 0(a 2 s ) correction to the A — >• gg decay width, extending the result of Refs. 0, [12| by 
one order. Our final result is listed in Eq. (p5|). For = 100 GeV, the overall QCD 
correction factor turned out to be as large as 1 + 0.68 + 0.23 = 1.91. It is slightly larger 
than the corresponding correction factor for the H —>■ gg decay width of the SM Higgs 
boson H with the same mass, which was found to be 1 + 0.66 + 0.21 = 1.87 [|Kj. The FAC 
PB| , PMS |3"4] , and BLM predictions for the 0(a 2 s ) correction to the A — > gg decay 
width are 0.38, 0.38, and 0.34, respectively. The corresponding results for the H — > gg 
case are 0.36, 0.36, and 0.33. 
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